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Abstract 

The mixing parameter ycp for the neutral D system can be obtained 
by comparing the lifetimes of CP eigenstate decays of the D mesons. 
Recent results on the lifetime differences from the Belle, CLEO, FO- 
CUS (E831) experiments are summarized in this article. The neutral D 
decay modes analyzed by the experiments are: The D° K--K+ mode 
with the assumption that it is an equal mixture of CP even and odd 
eigenstates, and the CP even modes Z)" — > and — * tt^tt^. 



1 Introduction 

Time dependent interactions of and decays can be coupled by the off- 
diagonal terms in the Schrodinger equation. If CP is conserved in neutral D 
decays, the eigenstates will be described as (Z?° ± I?o)/\/2 with mass and width 
as M ± 1/2AA/, T ^ 1/2 AT. Many calculations conclude that Standard Model 
contributions to x{= AAZ/F) and y(= Ar/2r) parameters are very small, and 
new physics effects can manifest themselves as an unexpected large measure- 
ment in the x parameter pQ. Experimentally, the mixing parameters x and y 
are accessible either by searching wrong sign decay of neutral D mesons or by 
comparing lifetimes of CP eigenstate decay channels. In this article, we summa- 
rize recent efforts of various experiments who are exploring the charm mixing 
sector using the second method. From now on, we assume that the charge 
conjugate decay modes give the same results as the original decay modes. 

For example, if we assume — > K~it'^ decay or its charge conjugate mode 
are advancing with equal amount of CP even and odd events, the width of the 
decay channels is described by T{Kti) (Fi -|-F2)/2. Meanwhile, D° K-K+ 
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and iD° ^ vr 7r+ decay modes are CP even and their width is described by 
Then ycp can be measured as: 

r{D ^ Kn) r{D ^ Kn) 

r^iJ KK) t(D — > tttt) 

In the following sessions, the analysis of the FOCUS experiment will be 
shown as well as the preliminary results from the Belle ^^'^ the CLEO 0] 
experiments. Implications of these measurements compared to the result of the 
wrong sign 13° decay analysis done by the CLEO experiment |5j will be discussed 
at the end. 



2 The FOCUS CP lifetime analysis 

The FOCUS (E831) collaboration consists of 100 physicists from USA, Italy, 
Brazil, Mexico and Korea. The experiment is a successor to the Fermilab 
E687 [H] experiment, and its main purpose is to study charm mesons and 
baryons produced by photon beams of about 180 GeV at a fixed target spec- 
trometer. During 1996-1997 run, the spectrometer collected more than 1 mil- 
lion reconstructed golden mode charm decays, D —^ Kn, D — > Knir, and 
D Ktttttt. The FOCUS group used two modes of neutral D decays, D Ktt 
and D —>■ KK, to obtain the ycp parameter. 

2.1 The FOCUS spectrometer performance 

The target material consists of 4 slabs of Be or BeO, depending on the data run 
periods. The segmented structure allows 62% of charm particles produced in 
the target material decaying in air, reducing the absorption correction necessary 
to the lifetime analysis to a small amount. Two sets of silicon microstrip vertex 
detectors are built to reconstruct vertices. The components of the first detector 
set interspace the target slabs, while the second set is located downstream of 
the target. The combination of the relativistic momentum of charged particles 
and excellent vertex detector performance produces an extremely good proper 
time resolution of 30 is for neutral D mesons, which corresponds to 8% of its 
lifetime. The good time resolution allows the FOCUS group to use a binned 
likelihood method free of deconvolution corrections to fit lifetime evolution of 
charm particles. 

The reconstruction algorithm of charm events is optimized to have a fiat 
acceptance over the fit variable, "reduced proper time" . The D production 
and decay vertices are reconstructed by a candidate driven algorithm. A cut 
is imposed on the amount of detachment (Z/cr) between the production and 
the decay vertices to suppress backgrounds coming from light quark events. 
Reduced proper time is a variable used by fixed target experiments to measure 
lifetime of particles, defined as the difference between the proper time and the 
amount of the detachment cut {t' = t — Na). 
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The FOCUS spectrometer has 3 stations of threshold Cerenkov detectors to 
identify charged particles. A continuous, likelihood-based identification system 
is used. This algorithm gives the flexibility useful in assessing systematic errors 
due to backgrounds. 

2.2 Event Selection and Fitting 

The event selection is optimized to obtain a sample which has a flat efficiency 
over the fit variable, t' . First, neutral D meson candidates are required to 
pass basic section cuts of detachments (Z/cr > 5) and for tracks reconstructed 
as kaons, the Cerenkov hypothesis to be a kaon should be favored over the 
hypothesis to be a pion by a factor of 7.39. Second, the candidate is either 
identified as a decay of D* (tagged sample) or required to pass more stringent 
particle identification cuts and kinematic cuts (inclusive sample). As a result, 
the FOCUS group obtained 119,738 D ^ Kir candidates and 10,331 D KK 
candidates. 

Histograms of 20 bins with 200 fs width are used to fit the reduced proper 
time of the D — + Ktt and D KK candidates. In each time bin, the amount 
of signal is estimated as those of exponential decays adjusted by the absorp- 
tion/acceptance correction function f{t') obtained from Monte Carlo (MC) sim- 
ulation. The amount of backgrounds is estimated from the sidebands of D mass 
lineshape. In case of the D KK sample, as shown in Figure ^ (b), a small 
amount of Z? — > Kn refiection background is still lingering in the sideband area. 
Since the time evolution of the reflection background in the D KK sam- 
ple is governed by the D — > Ktt lifetime, lifetime fitting should be done for 
both decay samples at the same time. The fit parameters are: D — > Ktt life- 
time {t{D Kit)), the mixing parameter ycp, and normalization parameters 
of background events for the D — > Ktt sample {Bk-k) and for the D KK 
sample {Bkk)- 

2.3 Results and Systematics 

Figure ^ (c) shows the reduced proper time evolution of the background sub- 
tracted D Ktt and D KK candidates in the final fit result. The lifetime 
and the lifetime differences of the neutral D mesons are fitted as follows: 



The systematic error in ycp is obtained by varying selection cuts and by chang- 
ing lifetime fit options. The systematic error on t{D — > Ktt) is not shown, 
since inputs from other neutral D decay channels, for example, D — > Ktttttt, are 
needed. 



ycp = (3.42 ± 1.39 ±0.74)% 
t{D ^ Ktt) ^ 409.2 ± 1.3 fs 



(2) 
(3) 
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3 The Belle lifetime study 



The Belle collaboration presented a preliminary measurements on lifetimes of 
and Ds mesons at the 4th International Workshop on B physics and 
CP Violation (BCP4) [2j. The neutral charm mixing study is shown as a part 
of their lifetime analysis package, based on the data set of ll/f& produced by 
an asymmetric e+e" collider, KEK-B, operating at 12 GeV T(4s). The mixing 
study is performed with two neutral D meson decay modes, — > K~'!r~^ and 
D° K-K+. 

3.1 Event Selection 

To suppress the bb background smearing into the charm sample, the momentum 
of D candidates is asked to be more than 2.5 GeV in the T(4s) frame. Iden- 
tification of charged particles is determined by using information from various 
detectors: a central drift chamber (CDC), a time of flight system (TOF), and an 
aerosol Cerenkov counter (ACC). K/n separation of up to 3.5 GeV is achieved. 
To enhance purity of the D sample further, a decay angle cut is applied for each 
mode. For signal events, the decay angle distribution [cosOo] should be flat. 
The requirement of the D* tag is optional. From the selection, the Belle group 
found 90, 601 ± 387 D Kit events and 7, 451 ±11S D ^ KK events. 

Proper lifetimes of charm mesons are calculated from the distance between 
the charm decay vertex and the production vertex {t — I ■ m{D) / p{D) ■ c). The 
decay vertex is reconstructed using all tracks from the charm candidate. The 
resulting D flight path is extrapolated to the beam interaction region and a 
production vertex is obtained. 

3.2 Lifetime and ycp fit results 

An unbinned maximum likelihood fit is used to obtain the charm lifetime and 
the mixing variable ycp- The likelihood function for each event is described 
by a sum of following terms: The signal term described by an exponential of 
lifetime tsig convoluted by a resolution function Rsig, a term for background 
events described by an exponential of lifetime tbg convoluted by a resolution 
function Rbg, and another term for background events which do not have a 
lifetime and convoluted by the Rbg- For signals (backgrounds), the resolution 
is 155 (160) is for 56 (60)% of events and 348 (350) fs for 44 (40)% of events. 

A combined lifetime fit using both D Kn and D KK distributions is 
performed and the following lifetime parameters are obtained: 



The errors shown above are statistical only. A detailed study on systematic er- 
rors is currently going on. The numbers are consistent with separately measured 
neutral D lifetimes; r(D Kn) = 414.5 ± 1.7 fs, t{D KK) = 409.8 ±6.3 fs 



t{D Kn) = 414.5 ± 1.7 fs 



(4) 
(5) 
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The Belle measurements on the lifetime and the mixing parameter ycp are 
consistent with the ones from other experiments. Figure |2] shows the mass and 
lifetime distributions used to extract the fit results. 



At the BCP4 conference, the CLEO collaboration summarized their recent ef- 
forts on measurements of D° — mixing, CP violation in £>° decays, and D*^ 
width 0]. They showed the preliminary analysis on lifetime using D — > Ktt, 
D KK, and D ^ tttt decay channels, based on a data set of 9.0/f6 produced 
by a symmetric e+e^ collider, CESR, operating at 10.6 GeV (T(4s)). 

4.1 Event Selection 

The kinematic characteristics of charm particles created in the CLEO detector 
are similar to those created in the Belle. The D* tag is required to select neutral 
D events, exploiting the excellent CLEO resolution on the mass difference be- 
tween D* and D" (ctq = 190 ± 2KeV). The CLEO group selected 20, 272 ± 178 
D Ktt, 2, 463 ± 65 D ^ KK and 930 ± 37 Z) tttt events. 

4.2 Lifetime and ycp fit results 

The CLEO analysis is also performed using unbinned maximum likelihood fits 
to proper time distributions of three D'^ decay channel candidates. The fitting 
method is optimized for finding lifetime differences, not for finding lifetime of 
neutral D mesons. The resolution function of proper time in the likelihood is 
described by three Gaussians: One Gaussian for events with correctly measured 
proper time and two Gaussians for events with mismeasured proper time. For 
most events, the resolution is about 35% of the neutral D lifetime. The following 
t{D^) values are obtained for each decay channel: r(Z?° K~7r~^) = 404.6 ± 
3.6 fs, t{D° K-K+) = 411 ± 12 fs, and t{D° tt-tt+) = 401 ± 17 fs. 
The t{D'^) values are highly correlated to the 1?° mass constrained in the fit 
procedure. This fit technique is not used in the dedicated D lifetime analysis 
published by the CLEO group 8 . Mixing parameter values are calculated as 
differences in lifetimes among three decay channels: 



This preliminary result is consistent with the measurements from the other ex- 
periments. Dominant systematic errors are coming from statistical uncertainty 
in MC lifetime study (9 fs), background description (8 fs), proper time resolu- 
tion model (5 fs), and fit procedure (5 fs). 



4 The CLEO L>" lifetime study 



2/CP (combined) = (-1.1 ± 2.5 ± 1.4)% 



ycp{KK) = (-1.9 ± 2.9 ± 1.6)% 
ycp{-KTT) = (0.5 ± 4.3 ± 1.8)% 



(6) 
(7) 

(8) 
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PDG2K 



412.6 ±2.8 11 



E791 


413 ± 3 ± 4 112 included in PDG2K 


CLEO 


408.5 ± 4.1 ± 3.5 8 included in PDG2K 


FOCUS (E831) 


409.2 ± 1.3 ±x* ^ 


SELEX (E781) 


407.9 ±6.0 ±4.3 [El 


Belle preliminary 


414.5± 1.7±a;* 



Average of 411.1 ±1.0 

recent values = 6-7 for dof = 4 



Table 1: The recent results on D° lifetime from various experiments are shown 
in chronological order. The average in the table is obtained from the most recent 
five measurements. The numbers are consistent between one another. The E791 
and CLEO values are already included in the PDG 2K edition JT]. x* denotes 
that corresponding systematic errors are not shown by the experiments yet. 

5 Conclusion 

The most recent five results on lifetime are shown on Tabled The mea- 
surements obtained by the FOCUS and the Belle experiments from the lifetime 
difference analysis are comparable to those from the dedicated D'^ lifetime anal- 
ysis. 

Table |31 shows compilation of ycp values obtained by various experiments. 
The numbers are comparable to one another, producing an average of (1.8 ± 
1.0)% with x^/dof = 2.3/3. On the other hands, the CLEO experiment 
published an impressive paper on the neutral D mixing parameters, x' and 
y' by analyzing lifetime evolution of wrong sign decays of neutral D mesons 
(D° — > K'^Tr~) 0. On the surface, the average ycp value is comparable to 
the CLEO wrong sign results, -5.8% < y' < 1% at 95%CL. But one has to 
be careful when combining ycp and y' together. The y' and ycp parameters 
are differed by a strong angle cos 5 originated from final state interactions. The 
strong angle has been considered to be small, but a recent theoretical paper 
suggested possibilities of not so small cos 5 [HI . The results on two y parameters 
have comparable errors but opposite signs, which prompted interesting discus- 
sions and speculations on neutral charm mixing |10| . This situation implies that 
mixing measurements in charm sector are challenging studies. We are eagerly 
waiting for up-to-date inputs from B-factories and next generation experiments 
to further enhance understanding in ever fascinating mixing phenomena. 
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E791 (0.8±2.9± 1.)% .12^ 

FOCUS (E831) (3.42 ± 1.39 ± 0.74)% 
Belle preliminary (1.16^111)% El 

CLEO preliminary (--1.1 ± 2.5 ± 1.4)% T 



Average of ycp (1.8 ± 1.0)% 

= 2.3 for dof = 3 



Table 2: The recent results on the lifetime difference ycp obtained by various 
experiments via D — > Kir^ D KK, and D tttt channels are shown in 
chronological order. The numbers are consistent between one another. One has 
to be careful if ycp values should be compared with y' parameter measurements. 
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Figure 1: (a) Reconstructed mass distribution of K^tt^ and its conjugate 

candidates from the FOCUS experiment. The yield is 119 738 K~Tr^ and K^n^ 
signal events; (b) That of D° K-K+ candidates. The yield is 10331 K~K+ 
signal events. The vertical and dashed lines indicate the signal and sideband 
regions used for the lifetime and ycp fits; (c) Signal versus reduced proper 
time for the D Ktt and D KK events shown above. Each data point is 
background subtracted and includes the (very small) Monte Carlo correction; 
The plots shown here are reprinted from P]. 
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Figure 2: Preliminary results from the Belle experiment on neutral D lifetime 
differences [3); (a) Reconstructed mass distribution of D° — > _?sr~7r+ and its 
conjugate candidates. The yield is 90, 601±387 K^-k^ and K^tt^ signal events; 
(b) The proper time distribution of the D — > Kn sample. The dashed line is 
the estimated background in the sample; (c) Reconstructed mass distribution 
of D° K-K+ candidates. The yield is 7, 451 ± 118 K- K+ signal events; (d) 
The proper time distribution of the D — > KK sample. The dashed line is the 
estimated background in the sample. 
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Figure 3: Preliminary results from the CLEO experiment on neutral D lifetime 
differences 0]; (a) Reconstructed mass distribution of K^tt^ and its 

conjugate candidates. The yield is 20, 272±178 K^tt^ and K^tt~ signal events; 
(b) The proper time distribution of the D Ktt sample; (c) Reconstructed 
mass distribution of Z?" R- K+ candidates. The yield is 2, 463 ± 65 R- K+ 
signal events; (d) The proper time distribution of the D — ^ KK sample; (e) 
Reconstructed mass distribution of 7r~7r+ candidates. The yield is 930 ± 

37 K^K^ signal events; (f) The proper time distribution of the D ^ tttt sample 
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